Interaction Between Kink and 
Radiation in 4 Model* 

Tomasz Romanczukiewicz^ 
<Vi ! Institute of Physics, 

q ! Jagiellonian University, Reymonta 4, Krakow, Poland 

(N : 

H ' February 1, 2008 

| Abstract 

> ' 

00 ■ The radiation from oscillating kink in (1+1) dimensional relativis- 

l/"") , tic cf) 4 model is considered. Both analytical and numerical approaches 

| are presented and the comparison between these methods is discussed. 

■ Acceleration of the kink in external radiation is calculated and numer- 

cn , ical results are also presented. 

jC ! 1 Introduction 

CD . 

^ ■ Kink is a simple example of topological defect. Topological defects i.e. 

>• . monopoles, vortices, domain walls play very important role in modern physics 

starting from flux tubes in QCD £Q, vortices in liquid helium and super- 
conductors, domain walls in magnetics, various defects in liquid crystals to 
strings in cosmology j2] (although recently it has been found cosmological 
strings are not as important as they had appeared to be .3j). Although 
these objects were studied for many years there are still many unanswered 
questions, mostly concerning their dynamics. There are some open problems 
about their interactions with other topological defects and external fields jlj . 
Even dynamics of their internal degrees of freedom is not fully understood 
and needs a lot of examinations In most cases we built our knowledge 

from experiments and numerical calculations. In theory, topological defects 
are described by nonlinear partial differential equations. Most of them one 
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can solve exactly only in special cases, for the rest of them one can apply 
only approximations. 

In this paper we focus on perturbation method around a well known ana- 
lytical static solution namely the kink. The results will be compared with 
numerical calculations. Very similar problems were discussed in [Zj, [Ej and 



Manton in [7j considered problem of kink and antikink creation and anni- 
hilation process and its relationship with the vibrational mode. He presented 
theory concerning the behavior of that mode and its decay. In this paper 
we develop methods presented in [Z| and apply to different initial conditions, 
more suitable for computer simulations. We verify Manton's and our pre- 
dictions via numerical calculations. In jS] there are another methods for 
calculating the radiation presented in a context of domain wall. The radia- 
tion from squeezed kink is investigated in jH]. 

Usually when one wants to find out whether the system is stable, one adds 
small perturbation and study its evolution. If the perturbation does not grow 
in time it means the solution is stable. Otherwise it is unstable. The problems 
of stability and relaxation process is well examined for dissipative systems 
(for example systems described by diffusion equations or wave equations with 
damping). The model considered in the present paper is energy conserving 
and therefore we can use stability only in a local sense. The relaxation 
process must be based upon radiation of the redundant energy into infinity. 
The static kink solution has the lowest possible energy in its topological sector 
(Bogomolny bound). This is a reason for us to treat kink as an attractor. 
The system with localized perturbation will tend to our kink solution in finite 
regions of the one-dimensional space. 

In the sections 4 and 5 of our paper we apply methods presented in the 
sections 2 and 3 to kink interacting with external radiation. Because the 
kink is transparent to radiation in linear order of approximation our usual 
intuition fails. For small amplitudes, the kink instead of being pushed away 
by radiation pressure is accelerating toward the source of radiation. 

In literature there are considered mostly kinks interacting with a constant 
external force [3] or a force oscillating in time ^U] or [TT] . 



Let us consider one dimension real scalar field theory described by the equa- 
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We use natural unit system (c = 1). There exist well known static solutions: 



4> s {x) = ±tanh(x — x a ) (2) 

which will be referred to as kink (with + sign) and antikink (with — ). With- 
out loosing generality one can choose x Q = due to translational invariance. 
Let us add a small perturbation to the kink solution: 

<) ) {t,x) = <f )s {x)+^t,x). (3) 

Eq. ((TJ) can be rewritten in the form: 

£ + L£ + iV(O=0, (4) 

where linear operator L has the form: 

L = -J^+[4-6(l-«], (5) 
and N denotes the part nonlinear in£: 

iV(0 = 6U 2 + 2£ 3 - (6) 

In the first approximation we assume that £ is small enough to neglect the 
term (JHJ). We seek solutions of the linearized eq. (jH) in the form £(t,x) = 
e luJt r](x). We obtain the eigenvalue problem (very similar to the Schrodinger 
equation) : 

Lr] = u 2r q. (7) 

The above equation has two solutions vanishing in infinity. One of them for 
uj = 

Vt(x) = — ^ (8) 
cosh x 

is called translational zero-mode because it is responsible for small transla- 
tions of the kink: <f) s (x + 5x) = <f) s (x) + 5xr] t (x) + 0(Sx 2 ). This mode plays 
a very important role when one considers the evolution of a system with an 
external force [1] or evolution of a system itself in more dimensions (domain 
wall) jHj or interaction between two or more kinks. The other normalizable 
solution is a vibrational mode with Ud = V^: 

Vd(x) = Sm ^ 2 X . (9) 
cosh x 

It is very important that in linear approximation this mode oscillates without 
loosing energy. It is a quasistationary solution because only couplings in 
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higher orders re responsible for decaying of this solution due to radiation. 
The energy is carried away by scattering modes in continuous spectrum: 

Vk ( x ) = e ikx (3 tanh 2 x - 3ik tanh x - 1 - k 2 ) , (10) 

where k is a wave vector: 

k 2 = u 2 -A. (11) 

We want to find the evolution of the system when initially only the discrete 
oscillating mode is excited. Let us construct a solution of the equation (jlj) 
substituting £(t, x) = A d (t)r] d (x) + rj(t,x) and assuming that A d and r\ are 
small we obtain the equation for 77 in 0(A 2 d ) order: 

^ + Lt7 + 6<M^ = 0. (12) 

This is an inhomogeneous linear equation for 77 with the source term T(t)g(x) 
where g(x) = r] d (x)(j) s (x) and T{t) = 6A d (t). Notice that it is proportional 
to the square of vibrational mode. We can make a time Fourier transform of 
the equation (fT2"J) and obtain: 

—cu 2 r](uj, x) + Lf/(u;, x) + T(cu)g(x) = 0. (13) 

It is an inhomogeneous linear equation for the spatial part but since we 
know the solutions ()10|) of the homogeneous equation we can easily construct 
appropriate Green's function in a standard manner: 

G k (x,y) = l~^ r]kiy)v - kix) X<V (14) 
(-wV-kiv^kix) x>y 

where W = —Hkik 2 + l)(k 2 + 4) is Wronskian of eq. (|13|). Therefore the 
solution has a form: 

T(uj) ( f x f°° 

fj(u,x) = — f r)- k (x) j dyr) k {y)g{y) + r) k {x) dyri^ k (y)g(y) 

'(15) 

We can use an asymptotic form of r]±k(x) ~ (2 — k 2 ± 3ik)e ±tkx for large x 
and obtain (the source is well localized around 0): 



After calculating this integral analytically we get: 

-/ x irk(k 2 + 4)(k 2 -2)f(uj) . . 

n(uj,x) = — r-rP- \ exp(-ikx) (17) 

/V ' 16sinh(f ) (2- k 2 -3ik) K ' K 1 

4 



OG 



Following Manton [7j we consider an oscillating mode A d {t) = acos(udt). 
Square of the amplitude is given by A? d = |a 2 (cos(2u; d t) + 1). Since constant 
solution carries no energy we take into account only the time dependent part: 
T[t) = 3a 2 e 2uJdt . Fourier transform has a form: T(u) = —6ira 2 5(uj — u ), 
where uj q = 2ua- Substituting the source term into fTTj) and calculating 
inverted Fourier transform we obtain: 

, , nk (k 2 + 4)(k 2 - 2) 9 ... , . , 

ri(t,x) = ; ° , ° " exp \i (u a t - k x . 18 

/V ; 32sinh(^) (2- k 2 -3ik ) PLV n K 1 

Since we are only interested in real part of the above equation the radiation 
has a form: 

r](t, x) = Qa 2 cos (cu t — k Q x + 5) , (19) 

where 

kl = - 4, 



and 5 is a phase which we are not interested in. The same radiation goes 
also to — oo. We have made numerical simulations for various a Q , where a Q 
is initial amplitude of the oscillating mode (in the following section we show 
that the amplitude decreases due to the radiation). In large distance from the 
kink we measure the outgoing field. The difference from the kink solution 
is sketched in figures ^ and 121 for a Q = 0.05 and 0.4 respectively (for the 
clarity we have plotted only the local extrema of the field instead of whole 
oscillations). The amplitude predicted by eq. (fT^j) seem to fit very well in 
the first figure. On the second Figure one can also see the damping which 
will be explained in the next section. 

As we can see from eq. |T9|) the outgoing radiation has only one single 
frequency which is twice the frequency of oscillating mode. This is true 
only in approximation, because the coupling in cubic term in (JH)) will result 
in frequency 3^- If one wants to test this procedure numerically one has 
to be very careful. The first idea is to pose initial conditions in the form 
0(0, x) = (p s (x) + arid(x) and 0(0, x) = 0. That means there is no radiation 
in t = 0, so there can be no source for t < 0. Source must have a form 
T{t) = 3a 2 Q(t)e luJot , where Q(x) is a Heaviside function. Fourier transform: 

f(u) =3a 2 ( -tt5{u - u ) + P — - — ). (21) 
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The field we calculate is equal to 



. . a* , r°° J irk(k 2 + A)(k 2 - 2) exp(i(ut - kx)) l 

■nit, x) = — P.V. / duj — J \ - ' — — — - -tt5(uj - uj ) + 

IK 1 2vr 32 sinh (f ) 2 - k 2 - 3ik) 



LU — LU r 
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Let us consider the second part of the integral in ()22|) containing factor 



I(t,x) = i / duF(u) 

UJ - UJr, 



(23) 



If k is a real number we have oscillations of the integrand along x axis (be- 
cause of the exp(— ikx) term), but when k is imaginary we have exponential 
decay, therefore for large x we need to take into account only u for which k 
is real, that is u > £1, where Q = 2. We can rewrite the above integral in a 
form: 

roc i(uj—uj )t 

I(t,x) =ie iuJat / duF(u) 

Jn 

After introducing v = [uj — u Q )t we obtain: 



UJ — UJ r 



(24) 



I(t,x) =%e Wot 
For large times we can expand F in |: 



dvF - + uj ) — 
(n-w„)t Kt J v 



(25) 



I(t,x) = ie iUot / dvF(uj ) — + F\uj )- 
Because Q — uj < and t tends to infinity the first term is equal to: 



(26) 



lim 

t— >oo 



dv F (uj ) — = F(uj ) / dv — = iirF(uj 



(Q — LV )t 



oo &lv 



(27) 



The second term we can integrate immediately after adding a term — ev to 
the exponent and obtain: 



I(t,x) = -e Wat 
But F(u) = e~ tkx f(uj) therefore 



TlFiuJn 



3 i(f2— uj )t 



-F'(uj n 



(28 



F'(u ) 



-ixf(uj ) 



dk 
duj 



+ f'M e 
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For large x we can neglect the second term: 



u e 



,i(Q—uj )t ' 



I(t,x) 



_ e iu a t 7T _ ix 




(29) 



t 



Finally the radiating part in our approximation has a form 



r](t, x) = Qa 2 1 — i 




exp [i(uj t — k Q x + §)] . 



(30) 



k Q lixtjx 



Apart from oscillation with frequency uj we have some decaying modulations 
with frequency uo D — Q m 1.46. Speed of the decay is determined by the point 
of observation x. The larger x the slower decay. Figure 3 shows the Fourier's 
transform of the numerically computed field in x = 30. As one can see there 
are three easily visible peaks. The largest one is at uo = uj = Iuj^. The 
next one is at another harmonic frequency 3ud- There is also one peak for 
uj = 1.15u)d = 2 = Q which is responsible for decaying modulations of the 
field. On the next Figure we have plotted envelope of oscillations evaluated 
from eq. (|22[). On Figure 5 one can see the behavior of the envelope for 
large times which seems to be well approximated by eq. (30). The measured 
frequency is really u — Q ~ 1.42. We also present the growth of the field 
for small times on figures 6 and 7. On the next pictures we have sketched 
numerically evaluated field (from eq. ()22j1 ) and, for comparison, numerically 
calculated solution of the full partial equation (1). Because of the numerical 
errors there are only the largest modulations visible. 

3 Backreaction 

In the above approximation we have found the field radiated to the infinity. 
We assumed that amplitude of the oscillating mode does not change in time. 
This is not true if the system is isolated. Radiation carries away energy 
and that leads to a decrease of the amplitude since the system is energy 
conserving. Let us calculate the rate at which the energy is being lost during 
radiation. For scalar field the energy escaping from a segment (xi, X2) is equal 
to ^ = <j>4>'\ xi ~ ^\ X2 - After taking 4>(t,x) = 4> s (x) + ar] d (x) cos(cu d t) + 
Qa 2 cos(u t — k Q x + 5) and averaging over a period we obtain the energy 
change inside large segment: 



dE 
~dt 




0.020a 4 



(31) 
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and since energy of the vibrational mode is E = a 2 we find a differential 
equation for the amplitude of that mode: 

d 4~ = -Q 2 k u a A . (32) 
at 

We can easily find a solution with an initial condition a(t = 0) = a Q : 

a{t) = . a ° (33) 
a/1 + Q 2 k u a 2 t 

From the above we obtain time dependence of radiation amplitude: 

C [x, a ) + Qk u t 

where 

C(x, a a ) = ——r - Qk u t (x) (35) 
Qa 2 D 

is constant in time. t a (x) is the time it takes for the field to travel from origin 
to point x where it is observed. 

This very precise prediction we can compare with purely numerical re- 
sults. Maxima and minima of the field with initial a Q = 0.4 measured for 
x = 200 are shown on Figure |2J There are also two hyperbolic function 
fitted to these data. The same procedure was applied for other amplitudes. 
The factor that stands besides t in eq. (|34jl in the theory is equal to 0.444. 
Numerical results show that its average value is 0.449 which is not signifi- 
cantly different. The value does not change much (within 2%) for different 
a Q (0.05^0.7). 

In FigureElwe have sketched fitted C versus a a . As one can see the fitted 
function in the form (|35jl is a very good approximation of the results. We 
find Q = 0.04607. In the theory above we have 0.0453 but it is still within 
numerical error. The last parameter in numerics is equal to t a = 336. The 
group velocity of the field is equal to v = §r|,_, = ^ ~ 0.816 and hence 
the distance x = 200 than theoretical value of t„jh = 245. In fact we actually 
can observe that the radiation cames to x — 200 in this time (figure Q and 
I2J). If one plots the fitted function and the theoretical one on one graph one 
can find they are almost indistinguishable, especially for small a Q where our 
theory is the most accurate. 



4 Anharmonic corrections 

As we saw in previous section, the field radiated from oscillating kink is of 
order O(A^) (eq. (fTHj) . Ad = acos(u i)). The oscillating mode reaction to 
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scattering modes will be of order 0(A d ) (in eq. (|3*3^1 a(t) ~ a — ^a Q 2 k ui t). 
But there is still nontrivial behavior of amplitude in the order 0(A 2 d ). Let 
us consider only the discrete mode. Substituting £(t, x) = A d (t)r] d (x) in eq. 
P|) in order 0{A 2 d ) we obtain the following equation: 

[A d + 3A d ) Vd + 6A> S ^ = Q (36) 

We are interested only in evolution of A d (t).We can rewrite the source term 
in a form 4> s f] d = a Vd + V±- To get rid of the orthogonal term t]± we project 
this equation onto r\ d . In order to do that we integrate both sides of the 
above equation with rj d and we obtain an ordinary differential second order 
equation (neglecting the perpendicular part): 

Ad + 3A d + %;A\ = 0. (37) 
lb 

We solve this using standard series method: A d = A^ + A^ + A^ + • ■ ■ . 
We take initial conditions: A(0) = 0, A(0) = a. In order O(A^) we obtain 

i(!) + 3A« = 0, (38) 
and the solution is A^ = a cos V3t. In second order we have: 

A (2) + 3A (2) + 97T fl 2 cog 2 ^ = q (3g) 

16 

The whole solution has a form: 

M*) = -% a2 + acos (^ t ) + cos(2V3t). (40) 

This function's maxima are equal to A max = a and minima v4 m j n = —a — 
||a 2 . Figure fTTl shows the oscillating mode evolution calculated numerically 
for a = 0.2. The field was measured in x = x max = arcosh \/2 ps 0.88, where 
the function r\ d has maximum equal \. The dashed lines represent \A max and 
\A m in calculated for that amplitude. As one can see our prediction works 
very well. 

Because the radiation depends on the square of the oscillating mode one can 
expect that minima are on more or less the same level while the maxima are 
changing (one is higher and the next one is lower and the following one is 
higher again). This effect is actually seen on Figure El 
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5 Kink in external radiation 



Let us now consider quite opposite process when kink is exposed to external 
radiation coming from +00: ae luJqt r] q (x) . Let us consider small amplitude 
limit. In the linear approximation there is no reflection in the potential of 
the soliton and therefore one cannot expect the system will behave as other 
similar systems in other branches of physics (there is no analogy to a particle 
exposed to an electromagnetic wave). In the first order the solution has a 
form (fTKJ) . where T(cu) = —3ira 2 5(uj — 2u q ) and g(x) = r] 2 (x)(j) s (x). The 
correction oscillates with frequency 2uo q . For x — > ±00 we can approximate 
the solution to form: 

fj(u,x^ 00) = — Iri-k(x) / dyri k {y)rj 2 q {y)(j) s {y) 

-T}- k {x) \ dyr] k (y)rj 1 q (y)(p s (y) + rj k (x) / dyri_ k (y)r) q (y)0 s (y) 



and 



f(u>) 



fj(w,x->-oo) = — (Vk(x) / dyr]^ k (y)rj q (y)(f) s (y) 

(42) 

-r} k {x) I dyr]^ k (y)rf(y)<p s {y) + r]_ k {x) I dyrj k (y)r] 2 (y)(j) s (y) 



where k = k(q) = (2uj q ) 2 — 4. In order to calculate the integrands with 
one finite but large limit we can use once again the asymptotic form of f] k {y). 
Integrands with both infinite limits can be calculated analytically via residua 
and finally we can write a simply expression for asymptotic behavior of the 
first correction: 

V+oo(t, %) = b + {u) q ) cos(u k t — kx + 5i) + c{u q ) cos(u k t + 2qt + S 2 ) (43) 

and 

?7_ 00 (t,x) = b-(uj q ) cos(u k t + kx — 5i) — c(u q ) cos(u k t + 2qt — 5 2 ). (44) 
where 

„ -60480 + 82224a;;? - 31620^ + 3707^ ± oA;(15120 - 11115a;? + 185 lui 

«.) = feir 



lOsinh ( — ^vr ) (Au 2 - 3)u 2 k 

(45) 
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cK) = ^(g 2 + iK 2 (46) 

and #1,2 are phases. The reflected part b+ for small frequencies is more or less 
the same order as 6_, but for large frequencies decays exponentially b + (u q ) ~ 
2 2 8 2 27 iiu q a 2 e~ 27TUJq . The transition part 6_ ~ ^u 2 a 2 . c is an amplitude of a 
wave coming from oo and going to — oo. This wave has exactly the same 
form at both sides and hence passes no energy nor momentum to the kink 
as one can see in details in the following calculations. 

For small x we do not know an analytical solution. The solution (which 
for small x can differ from rj± k ) oscillates around the kink solution with a 
frequency of the source wave 2io q . Because of the energy gathered in these 
oscillations kink as seen from a distance gains some extra mass. In following 
calculations we will be refering to that effect mass as M*. The bare mass of 
the soliton equals M = |. As mentioned in a section 3 the change of energy 
inside a segment (sci,^) equals ^ = </></>' |^ 2 . In our case we have from the 
right side <p(x) = <fi s (x) + A cos(u q t + qx) + b + cos{u k t — kx) + ccos(u k t + 2qt). 
The energy flowing into the segment averaged over a period is 

dE r 1 o 7 7 2 2 

— = -quJ q A - -ku k b + + qu k c . 

From the left side the field is <f>(x) = 4> s (x)+B cos(ui q t+qx)+b^ cos (2u q t + kx)- 
ccos(2u k t + 2qt). The energy is equal than: 

— = ~-quJ q B - -kuj k b_ - quj k c . 

The energy conservation takes a form: 

\qu q A 2 - l -qu q B 2 - \ku k {bl + b\) = M*^. (47) 

Where M*j = is an energy of the soliton. Apart from M* which is 

probably different from M we have two unknown kinetic variables B and 
7. The second one is responsible for the acceleration of the kink due to 
interaction with the radiation. We need one more equation. We can use the 
momentum conservation law: 

dPr 1 2 a 2 1 1 2j2 . 1 2 2 

ir = r« A + r kh + + r kC 

and from the left side 

dPl 1 2 r>2 1 2j2 1 2 2 

— — = — oj n B ujJ) u h c 

dt 2 q 2 k ~ 2 k 
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and hence the second equation 



~ \^B 2 - y (bl - b\) = M*|M = F* (48) 

These equations are correct only for v = because the moving kink 
experiences the Doppler's effect and there should be corrections in u q ^, k 
and q. But substituting v = to above equations we obtain the solutions 
very easily: 

' B 2 = A 2 - 2- (b 2 _ + bl ) 
q v ' 

(49) 

uj 2 

F* = — ({k - 2q)b 2 _ + (k + 2q)bV) 
\ q 

We used the relation Uk = 2u q . The factor before 6_ is 

k-2q = 2(y/u*-l-y/u*-4) 

is always positive and therefore the soliton accelerates in the direction from 
which radiation comes. The radiation presure is negative! The force desribed 
by above equation is shown on Figure fT2l 

Because b + (oj q ) is usually very small, for certain uj q when b_(u q ) = the ac- 
celeration derived in this order of approximation vanishes and nonlinearities 
of higher orders play the crucial role. Numerical results show the soliton will 
be pushed by the radiation. The conclusion is that the soliton in most cases 
will stem the current. If the soliton reaches the speed when the frequency of 
the source wave due to the Doppler's effect uj' = 7 (oj q + v^Ju^ — 4) is the 
frequency for which acceleration vanishes the soliton moves with more or less 
constant velocity. Figure El shows the paths of zeros of the field (kinks) for 
frequencies around the minimum of the pulling force. If the source wave has 
a frequency very close to the minimum of the force the kink's acceleration is 
small and the soliton moves with the current. 

If the source wave amplitude is large the above approximation fails. There 
exists certain critical amplitude, depending on a frequency, for which the 
soliton is pushed away by the incoming wave (Figure EJ). 
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6 Conclusions 



In this paper we have investigated the behavior of excited oscillating mode 
and its decay using both analytical and numerical methods. We have proved 
conformity between these two methods. We have studied the interaction of 
the kink with external radiation field. We have found that in most cases the 
kink will be pulled by the radiation. We have explained the peculiar behavior 
of this system. Finally we presented numerical results for large amplitudes 
of the source wave when the kink is pushed away by the wave. 
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Figure 3: Fourier transform of the radiation (a = 0.8) 




500 



Figure 4: Envelope of the radiation 
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Figure 6: Envelope for small times 
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Figure 7: Envelope for very small times 
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Figure 8: Radiation calculated from integral (7) 
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Figure 9: Radiation computed from the full partial equation (1) 
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Figure 10: C(x,a ) defined by eq. (}35|) vs initial amplitude of the discrete 
mode; numerical data and fitted function. 
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Figure 11: Anharmonical oscillation of the discrete mode. Dashed lines are 
calculated theoretically. 
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Figure 12: Pulling force of the radiation vs frequency 
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Figure 13: Zeros of the field for different frequencies of the source wave and 
for the amplitude a = 0.2 
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Figure 14: Zeros of the field for uo = 3.5 and for different amplitudes 
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